Search for: All records

The paper is concerned with efficient time discretization methods based on exponential integrators for scalar hyperbolic conservation laws. The model problem is first discretized in space by the discontinuous Galerkin method, resulting in a system of nonlinear ordinary differential equations. To solve such a system, exponential time differencing of order 2 (ETDRK2) is employed with Jacobian linearization at each time step. The scheme is fully explicit and relies on the computation of matrix exponential vector products. To accelerate such computation, we further construct a noniterative, nonoverlapping domain decomposition algorithm, namely localized ETDRK2, which loosely decouples the system at each time step via suitable interface conditions. Temporal error analysis of the proposed global and localized ETDRK2 schemes is rigorously proved; moreover, the schemes are shown to be conservative under periodic boundary conditions. Numerical results for the Burgers' equation in one and two dimensions (with moving shocks) are presented to verify the theoretical results and illustrate the performance of the global and localized ETDRK2 methods where large time step sizes can be used without affecting numerical stability. 

In this work, we present a numerical method that provides accurate real-time detection for the widthsof the fractures in a fractured porous medium based on observational data on porous medium fluidmass and velocity. To achieve this task, an inverse problem is formulated by first constructing aforward formulation based on the reduced fracture model of the diffusion equation. A parameter estimationproblem is then performed online by utilizing a direct filter method. Numerical experimentsare carried out to demonstrate the accuracy of our method in approximating the target parameters. 

In this paper, we present efficient numerical schemes based on the Lagrange multiplier approach for the Navier-Stokes equations. By introducing a dynamic equation (involving the kinetic energy, the Lagrange multiplier, and a regularization parameter), we form a new system which incorporates the energy evolution process but is still equivalent to the original equations. Such nonlinear system is then discretized in time based on the backward differentiation formulas, resulting in a dynamically regularized Lagrange multiplier (DRLM) method. First- and second-order DRLM schemes are derived and shown to be unconditionally energy stable with respect to the original variables. The proposed schemes require only the solutions of two linear Stokes systems and a scalar quadratic equation at each time step. Moreover, with the introduction of the regularization parameter, the Lagrange multiplier can be uniquely determined from the quadratic equation, even with large time step sizes, without affecting accuracy and stability of the numerical solutions. Fully discrete energy stability is also proved with the Marker-and-Cell (MAC) discretization in space. Various numerical experiments in two and three dimensions verify the convergence and energy dissipation as well as demonstrate the accuracy and robustness of the proposed DRLM schemes. 

This paper is concerned with the numerical solution of the flow problem in a fractured porous medium where the fracture is treated as a lower dimensional object embedded in the rock matrix. We consider a space-time mixed variational formulation of such a reduced fracture model with mixed finite element approximations in space and discontinuous Galerkin discretization in time. Different spatial and temporal grids are used in the subdomains and in the fracture to adapt to the heterogeneity of the problem. Analysis of the numerical scheme, including well-posedness of the discrete problem, stability and a priori error estimates, is presented. Using substructuring techniques, the coupled subdomain and fracture system is reduced to a space-time interface problem which is solved iteratively by GMRES. Each GMRES iteration involves solution of time-dependent problems in the subdomains using the method of lines with local spatial and temporal discretizations. The convergence of GMRES is proved by using the field-of-values analysis and the properties of the discrete space-time interface operator. Numerical experiments are carried out to illustrate the performance of the proposed iterative algorithm and the accuracy of the numerical solution. 

In contrast to the classical Allen-Cahn equation, the conservative Allen-Cahn equation with a nonlocal Lagrange multiplier not only satisfies the maximum bound principle (MBP) and energy dissipation law but also ensures mass conservation. Many existing schemes often fail to preserve all these properties at the discrete level or require high regularity in time on the exact solution for convergence analysis. In this paper, we construct a new class of low regularity integrators (LRIs) for time discretization of the conservative Allen-Cahn equation by repeatedly using Duhamel's formula. The proposed first- and second-order LRI schemes are shown to conserve mass unconditionally and satisfy the MBP under some time step size constraints. Temporal error estimates for these schemes are derived under a low regularity assumption that the exact solution is only Lipschitz continuous in time, followed by a rigorous proof for energy stability of the corresponding time-discrete solutions. Various numerical experiments and comparisons in two and three dimensions are presented to verify the theoretical results and illustrate the performance of the LRI schemes, especially when the interfacial parameter approaches zero. 

ABSTRACT This paper is concerned with efficient and accurate numerical schemes for the Cahn‐Hilliard‐Navier‐Stokes phase field model of binary immiscible fluids. By introducing two Lagrange multipliers for each of the Cahn‐Hilliard and Navier‐Stokes parts, we reformulate the original model problem into an equivalent system that incorporates the energy evolution process. Such a nonlinear, coupled system is then discretized in time using first‐ and second‐order backward differentiation formulas, in which all nonlinear terms are treated explicitly and no extra stabilization term is imposed. The proposed dynamically regularized Lagrange multiplier (DRLM) schemes are mass‐conserving and unconditionally energy‐stable with respect to the original variables. In addition, the schemes are fully decoupled: Each time step involves solving two biharmonic‐type equations and two generalized linear Stokes systems, together with two nonlinear algebraic equations for the Lagrange multipliers. A key feature of the DRLM schemes is the introduction of the regularization parameters which ensure the unique determination of the Lagrange multipliers and mitigate the time step size constraint without affecting the accuracy of the numerical solution, especially when the interfacial width is small. Various numerical experiments are presented to illustrate the accuracy and robustness of the proposed DRLM schemes in terms of convergence, mass conservation, and energy stability. 

The objective of this paper is to develop efficient numerical algorithms for the linear advection-diffusion equation in fractured porous media. A reduced fracture model is considered where the fractures are treated as interfaces between subdomains and the interactions between the fractures and the surrounding porous medium are taken into account. The model is discretized by a backward Euler upwind-mixed hybrid finite element method in which the flux variable represents both the advective and diffusive fluxes. The existence, uniqueness, as well as optimal error estimates in both space and time for the fully discrete coupled problem are established. Moreover, to facilitate different time steps in the fracture-interface and the subdomains, global-in-time, nonoverlapping domain decomposition is utilized to derive two implicit iterative solvers for the discrete problem. The first method is based on the time-dependent Steklov–Poincaré operator, while the second one employs the optimized Schwarz waveform relaxation (OSWR) approach with Ventcel-Robin transmission conditions. A discrete space-time interface system is formulated for each method and is solved iteratively with possibly variable time step sizes. The convergence of the OSWR-based method with conforming time grids is also proved. Finally, numerical results in two dimensions are presented to verify the optimal order of convergence of the monolithic solver and to illustrate the performance of the two decoupled schemes with local time-stepping on problems of high Péclet numbers. 

Abstract The Allen‐Cahn equation satisfies the maximum bound principle, that is, its solution is uniformly bounded for all time by a positive constant under appropriate initial and/or boundary conditions. It has been shown recently that the time‐discrete solutions produced by low regularity integrators (LRIs) are likewise bounded in the infinity norm; however, the corresponding fully discrete error analysis is still lacking. This work is concerned with convergence analysis of the fully discrete numerical solutions to the Allen‐Cahn equation obtained based on two first‐order LRIs in time and the central finite difference method in space. By utilizing some fundamental properties of the fully discrete system and the Duhamel's principle, we prove optimal error estimates of the numerical solutions in time and space while the exact solution is only assumed to be continuous in time. Numerical results are presented to confirm such error estimates and show that the solution obtained by the proposed LRI schemes is more accurate than the classical exponential time differencing (ETD) scheme of the same order. 

This paper is concerned with the numerical solution of compressible fluid flow in a fractured porous medium. The fracture represents a fast pathway (i.e., with high permeability) and is modeled as a hypersurface embedded in the porous medium. We aim to develop fast-convergent and accurate global-in-time domain decomposition (DD) methods for such a reduced fracture model, in which smaller time step sizes in the fracture can be coupled with larger time step sizes in the subdomains. Using the pressure continuity equation and the tangential PDEs in the fracture-interface as transmission conditions, three different DD formulations are derived; each method leads to a space-time interface problem which is solved iteratively and globally in time. Efficient preconditioners are designed to accelerate the convergence of the iterative methods while preserving the accuracy in time with nonconforming grids. Numerical results for two-dimensional problems with non-immersed and partially immersed fractures are presented to show the improved performance of the proposed methods.